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Abstract. In a cosmological context dust has been always poorly understood. That is 
true also for the statistic of Gamma-Ray Bursts (GRBs) so that we started a program 
to understand its role both in relation to GRBs and in function of z. 

This paper presents a composite model in this direction. The model considers a rather 
generic distribution of dust in a spiral galaxy and considers the effect of changing some 
of the parameters characterizing the dust grains, size in particular. We first simulated 
500 GRBs distributed as the host galaxy mass distribution, using as model the Milky 
Way. If we consider dust with the same properties as that we observe in the Milky Way, 
we find that due to absorption we miss ^ 10% of the afterglows assuming we observe the 
event within about 1 hour or even within 100s. 

In our second set of simulations we placed GRBs randomly inside giants molecular 
clouds, considering different kinds of dust inside and outside the host cloud and the 
effect of dust sublimation caused by the GRB inside the clouds. In this case absorption 
is mainly due to the host cloud and the physical properties of dust play a strong role. 
Computations from this model agree with the hypothesis of host galaxies with extinction 
curve similar to that of the Small Magellanic Cloud, whereas the host cloud could be 
also characterized by dust with larger grains. Unfortunately, the present statistics lack 
solid grounds, being based on hardly compatible observations, at different time from the 
burst and with different limiting magnitudes. To confirm our findings we need a set of 
homogeneous infrared observations. The use of coming dedicated infrared telescopes, like 
REM, will provide a wealth of cases of new afterglow observations. 
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Table 1. Physical characteristics of H2 galactic clouds. 
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1. Introduction 

The observations show that about 50% of the detected GRBs are not visible at optical 
wavelengths. Statistics refers to well localized bursts (Bloom, Kulkarni & Djorgovski 
2002) with an X ray flux rather similar to those that have been detected also in the 
optical band. The rapid decline of the optical afterglow (Fynbo et al. 2002, Berger et al. 
2002) can not explain the observations as well and that is why we call them dark bursts 
and an explanation has yet to be found. 

In two cases, GRB 970828 and GRB 990506, in spite of lacking the optical afterglow 
we were able to identify the host galaxies (z=0.958 and z= 1.3 respectively) so that in 
these cases it is evident that high redshift is not the cause of the optical flux extinction. 
A possible explanation is that the events have been obscured by dust. After the first 
part of our work was completed (Vergani 2002) Reichart (2001) and Rcichart & Price 
(2002) made a case for a strong absorption occurring in the molecular clouds where 
the event originates. This is in agreement with the current thinking that GRBs are 
related to massive star formation which is strongly correlated with very dusty regions. 
Whether we are dealing with a particular type of galaxies is not known. Ramirez-Ruiz, 
Threntham & Blain 2002 assert that we might be dealing with ULIRG (Ultra Lmninous 
Infrared Galaxies) or alike objects and that the extinction is essentially due to the dust 
distribution present in these galaxies. We do not have observational evidence that this 
is the case however. No matter how the dust is at work, we must also account for the 
local process of dust sublimation and understand how much dust a burst is capable of 
sublimating and sweeping out. 

We decided to tackle the problem first theoretically, and this paper report part of 
our work in this direction, and observationally by building the robotic NIR Telescope 
REM (Zerbi et al. 2002, Chincarini et al. 2003). The statistics of GRBs wiU also largely 
increase as soon as the Swift satellite will be launched (Gehrcls et al. 2003). 

In section 2 wc discuss how we model the dust itself, Section 3 is devoted to the 
construction of the simple host galaxy; a basic GRB distribution models and results are 
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presented in section 4. In section 5 and 6 we associate GRBs with giant molecular clouds 
and we explore its implications. Conclusions are drawn in section 7. 



2. Dust model 

The afterglow radiation reaches the observer after interacting with circunburst material, 
host ISM, IGM and the ISM of our galaxy. In this work we consider only the interaction 
with the dust in the ISM of the host galaxy. 

Our dust model is based on Mathis, Rumpl & Nordsiek (1977) model, improved by 
Mathis (1986). We suppose that dust is made by spherical grains composed by 50% of 
graphite and 50% of silicates with a grain size a distributed as 

n{a) = noa~^-^ (1) 

with amin — 0.005/im, Umax = 0.25/um and no, proportional to the neutral hydrogen 
density (n^f [cm^"^]), used bv lVenemans fc Blainl 2001. 

The variation in magnitudes produced on flux after interacting with dust grains is 

/>oo 

Am oc / Qej;tna^n{a)da (2) 
Jo 

with Qext ^extinction efficiency (sum of scattering and absorption efficiency, see 
I Van de Hulst 1957|l. 

By varying in our model the fraction of the two types of materials, the exponent of the 
grain size distribution and the grain size, the result is a change in the scale of the curve 
without modifying its shape. We cannot therefore disentangle their contribution from 
different amounts of hydrogen column densities (NH[cm^^]) along the line of sight. On 
the contrary the size of the grains play a strong role. Large grains, which could be present 
in regions with intense star formation ( |Maiolino, Marconi fc Oliva 2001| ) or in circunburst 
environment after the burst has sublimated the smaller grains IjVenemans fc Blain|l . cause 
the extinction curve to flatten. For dust composed only by large grains {amin = l/^m, 
o-max = 2fim) our computations predict a fixed value of Am ~ 0.1, in agreement with 
the theory, that predict for the case 2a > A a fixed value of Qext ^ 2 (Fig.^. The effects 
of modifying amin with fixed amax are not relevant. 



3. Galactic-like model 

We model an host galaxy similar to our Milky Way, where the dust is distributed as the 
neutral hydrogen, both molecular (II2) and atomic (HI). 

For the galactic hydrogen surface density we adopt the values taken from 



Scoville (1992), and Binney fc Merrifield (1998) 
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Fig. 1. Extinction curves produced by varying the dust grain size in case of an hydrogen 
cohimn density of lO^^cm^^ . Dotted steeper curve represents the extinction considering 
our Galactic-hke dust model (STD, see Tab.^ based on Mathis (1986). Dashed and sohd 
hnes describe the extinction curve in case of a dust composed by larger grains (OLG and 
LRG, respectively). 



For H2 , two gaussian fits were drawn on nuclear and disk data: namely the gas density 
has been expressed as: 



NH(bulge) oc e^I^^cm""^ (3) 
and 



-(H-5.03) 
— -2 



NH(disk) oc e ^d.s)^ cm"^ (4) 

for the bulge and the disk, where R is the radial distance in cylindrical coordinates (kpc) . 
The total mass of galactic H2 is distributed with 4 • lO*Af0 in the bulge and 2 ■ lO'^M© 
in the disk. The H2 clouds are present in two distinct morphologies, which we catalogue 
as giant molecular clouds (GMC) and dense clouds (DC) whose adopted characteristics 
are summarized in Tabled An algorithm assigns positions, with a random distribution 
weighted by Galaxy mass density, to GMC and DC. 

Neutral atomic hydrogen is present in a diffuse form throughout the disk system 



with a total mass of 4.3 • 10 0. Following Binney & Merrifield (1998) we simplify its 



distribution as being sharply confined in the radial interval 3 < i? < 18 (i? in kpc) 
with a spatial density represented by ni/[cm^^] = 0.797e^'* ^'^■'^^ (where h is the height 
from Galactic plane in kpc) . In Fig. [21 we plot the adopted fits for HI (dashed line) and 
II2 (solid lines) surface densities, while Fig. 13 shows the appearance of the host galaxy 
hydrogen distribution . 
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Fig. 2. Milky Way hydrogen column density as a function of radial distance R in kpc. The 



two solid curves are the two gaussian fits to the distribution reported in Scoville (1992), 
for the molecular hydrogen (H2 ) . Dashed line is the atomic hydrogen (HI) surface density, 
simplified as a box function between 3 and 18 kpc from galactic center. 



4. Basic GRBs absorption model 

We suppose that the GRBs are distributed as the luminosity (barionic mass) 
of the host Milky Way-like galaxy. Using the Milky Way photometric model of 



Kent, Dame & Fazio (1991) we have 



j,(i?,/,) = -^e-«/«-N/''°, 



with 



0.165 + 0.21(i?/i?o - -) 



kpc. 



(5) 



(6) 



Rd =3kpc, Rq =8kpc and the central surface brightness Id =1000 Lqpc ^. 
The bulge has been modeled using a King profile 

/(r)oc[l + (r/re)2]-3/2 



(7) 



with Tc =0.25kpc. 

The fraction of GRBs located in the Galactic disk and in the bulge is derived assuming 
a mass ratio disk to bulge of 4.5. 

The burst occurs randomly using as weighting function the distribution of mass of 
the galactic model. The observer is located randomly over a An solid angle. 
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Fig. 3. A rendered image of our host galaxy model. The top figure is the complete model 
with the sphere on the right which represents a randomly positioned GRB. Its line of sight 
crosses the galactic plane. Below the distribution of the atomic hydrogen (HI) and H2 
clouds (both GMC and DC) are separately shown. In the bottom picture the simulated 
distribution of the 500 GRBs. 

Each molecular hydrogen cloud is accounted for absorption and the amount of NH is 
integrated along the line of sight. The quantized NH amount is then 4.5 • 10^'^cm^'^ for 
the GMC and 2.25 • Hfi^croT^ for the DC. 

The atomic diffuse hydrogen is summed integrating along the line of sight from the 
GRB to a radius of 20 kpc from the galactic center. The HI contribution ranges therefore 
from a null value to 7 • lO^^cm"^, with a pick at about 5 • lO^^cm^^ not sufficient to 
obscure the afterglow. On the other hand the encounter of a single cloud yields values 
of NH of about 5 • lO^^cm"^ (GMC) or 5 • lO^^cm^^ (qq^ largely enough to completely 
absorb the optical afterglow. 

The total amount of NH is computed and catalogued for the whole 500 GRBs set 
and for the subset of 91 nuclear bursts. In Fig.0|the histogram of the total NH distribu- 
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Fig. 4. Hydrogen column density distribution for 500 simulated GRBs. 

Table 2. Extinction limits computed on the basis of the limiting magnitudes of different 
instruments and of GRB990123 afterglow light curve. NH values are computed using the 
dust model in Section 2. 
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(3) (4) In these columns are reported light curve values GRB990123-like. 

(5) (6) Prompt observation simulated with REM | |Vergani 2002| Table 2.3), late observation 
simulated with ISAAC and FORS (see |http: / /www.eso.org/oserving/etc/ 1 for its ETC time 
estimates) 

(9)(10) in units of lO^^cm-^ 



tion is shown for the two populations. To better appreciate the difference a cumulative 
distribution for a smaller range of column densities is reported aside. 

The amount of dark GRBs due to dust absorption can now be estimated, once we 
have fixed the limiting magnitude of our telescopes in the various passbands and the 
typical apparent magnitude of GRB afterglows. We associate randomly to each GRB a 
jet opening angle 6 following the law oc Q^'^-^^ that we have extrapolated from the known 
jet angles reported by Bloom, Frail & Kulkarni (2003) (see Fig[71 upper panel). Note 
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Table 3. Optical (R) and infrared (K) transient loss (%). Prompt and late observation 
cases are reported for z—1 galaxy position. 
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that the 9 distribution is truncated as 0.05 < 6[rad] < 0.6. Considering the luminosity L 
proportional to 9^^, we calculate the magnitudes R and K of GRBs at 100s and 5000s, 
taking as reference GRB990123 shifted at z=l (that is the observed mean redshift of 



GRBs, [Hurley et al. 2002| ), the color data by [Simon et al. (2001)|and using the relation 



mcRB "1990123 + 5 Iog(6'/6'990123) (8) 

where 6'99oi23 and 771990123 are the jet angle and the magnitudes of GRB990123. To 
compute the limiting magnitudes we assume to use REM telescope for the fast response 
to the GRB alert in R and K bands, FORS and ISAAC camera of ESO/VLT for the long 
term observations (Tab.[2J). 

For each GRB we calculate the extinction in R and K caused by the traversed NH 
column density. We are then able to compute the percentage of lost afterglows due to 
dust absorption under the assumption that the host galaxy at z—1 has dust properties 
similar to Milky Way. 

In Table|2lR and K magnitudes of prompt and late observations for GRB990123-like 
events are reported. The values of NH quoted in columns (9) and (10) represent the 
computed amount of NH needed, according to our dust model, to obtain the extinction 
values of column (7) and (8). The summary of Table |21 is a clear indication that only a 
relatively small fraction of GRBs is not observed due to dust extinction in the case of a 
host galaxy at z=l. 



5. GRBs in molecular clouds 

Our results show that, if the GRBs distribution is similar to the mass distribution and if 
the dust of the host galaxies has the same characteristic of Galactic dust, the percentage 
of dark GRBs due to dust obscuring is rather low. There is no relevant difference between 
R and K observations and between percentages relatives to observations taken at 100s 
and 5000s. It also happens that no burst, out of the 500 considered, occurs statistically 
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Fig. 5. Geometry around the GRB location explaining the shapes of the four regions 
whose content in terms of dust is parameterized in our models. The regions are spherical 
shells centered on the GRB location delimited by (the radius up to which all dust grains 
are destroyed), i?s(the radius up to which grains with radius smaller 1/im are destroyed), 
Rc (the radius of the giant molecular cloud) and the border of the galaxy. 6 is the jet 
opening angle associated to the GRB event. 

inside a molecular cloud. In this scenario the majority of dark bursts could be due to 
high redshift Ly-a absorption. 

We then consider the case that GRBs follow the distribution of giant molecular clouds. 
Physically the assumption is that of a strong connection between GRBs and massive stars 
formation. 

Within this framework, we place 5000 GRBs inside our modeled giant molecular 
clouds, that have similar characteristics respect the supposed typical GMC GRBs host 



described by Galama & Wijers (2001) 



The first aspect to consider is dust sublimation by the optical-UV flash accompanying 
the GRB, phenomenon likely confirmed by the observations of the afterglow light curve of 



GRB 990123 l|Akerlof et al. 1999|l . We consider the results of |Waxman fc Draine (2000) 
and lReichart 2001bl that compute the radius up to which the dust is sublimated, which 
is also a function of dust grain size. In case of a canonical distribution of graphite 

1/2 

and silicate grain size, the sublimation radius is Rs — 10^43 where ^49 is the 1- 
7.5eV (I6OO-I2OOOA) isotropic-equivalent peak luminosity of the optical flash in unit of 
IQ'^'^ergs^^ , that is the 1-7. 5eV isotropic-equivalent peak luminosity of the optical flash 
of GRB990123. 

Inside the host cloud we consider both the case of a standard galactic dust and the 
case of a dust, already present before the sublimation took place, with larger grain size 
with amin = 0.005/^m, amax = 2/xm, because the connection of GRBs and intense star 
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formation regions (Galama & Wijers (2001) I, as described in §2. Moreover, in the case 
of large grains, we have to consider two different sublimation radius: an inner one (r^) 
up to which all grain are destroyed and a larger one (Rs) up to which only the grains 
whit radius smaller than 1/im can be sublimated IjVenemans fc Blainjl . 

Rs and rg vary with the intensity of the peak luminosity of the optical flash, that we 
can suppose depending on the GRB jet opening angle 9. We calculate for each GRB Rs 
and Ts as Rs = i?s990i23 * ^'990123/6' and = rs99oi23 * 6*990123/6', where i?s990i23, '^3990123 
and 6*990123 are the two sublimation radius and the jet opening angle of GRB990123. 

To each of the 5000 simulated GRBs is associated a random line of sight passing 
through 4 regions: 9?1 is the inner region from the GRB to (the radius up to which all 
dust grains are destroyed); 5R2 goes from rg to Rs, where large grains (radii from 1/im to 
2 ^m) are allowed to survive; 5R3 is the rest on undisturbed host cloud dust; 5R4 is the 
host galaxy outside the host cloud, consisting of all other molecular clouds and of the 
diffuse medium (see Fig|SJ. In these regions we place 5 kind of dust yielding ten different 
models. 

In the first three models, to calculate the extinction in the region outside the host 
molecular clouds, we use our Galactic-like dust model. In 5R3, that is the part inside the 
host molecular clouds where there is no sublimation, 2 scenarios are explored: standard 
galactic dust (model 2 and 3) or a similar one but with larger grains {amax — 2/im, 
model 1). Moreover, in model 2 we suppose that dust in the whole cloud is standard and 
thus sublimated both in 5R1 and 5R2, whereas in model 1 and 3 the circunburst medium 
is supposed to be characterized by dust with Umax — 2/Ltm, thus in 5R2 dust grains larger 
than 1/zm are left. 

Hyorth et al., 2003 have recently supported the hypothesis that GRBs environment 
has an extinction curve similar to the one of Small Magellanic Cloud (SMC), so, in model 
Ismc, 2smc and 3smc we present the same scenario of model 1, 2 and n 3 but considering 
SMC extinction inside the cloud and/or in the host galaxy. Both the cases of extinction 
caused by larger grains or SMC-like dust agree with the non detection of the A — 2 175 A 



bump ( Galama & Wijers (2001) I, that is instead a significant feature of Galactic-like 
dust. 

Simulations by Perna et al. 2003, show that if we consider dust with properties similar 
to the Galactic one, burst energy sublimes very quickly silicate grains present in the 
circunburst medium differently from graphite made grain that need a longer time to be 
destroyed and whose sublimation is more efficient on the smaller grains. On this basis, 
we think reasonable to consider models (model lb, 3b, Ismcb and 3smcb) in which dust 
in 5R2 region is composed only by graphite grain with amin — 1/^m, Umax = 2/im. Table 
01 shows all the cases. 



Vergani, Molinari, Zerbi, Chincarini: Dust and dark GRBs 11 

Table 4. Regions around GRB location and their contents in term of dust type as adopted 
in our computations. Dust model used are: STD (our model of standard Galactic dust, 
with grain radii from 0.005/im to 0.25/im); LRG (with added larger grains, from 0.005/im 
to 2/im); OLG (only large grains, left from partial sublimation, with radii from 1/im to 
2/im); GRA (only large graphite grains, left from partial sublimation) SMC (dust with 
low extinction similar to Small Magellanic Clouds); SUB (no dust, completely sublimated 
by the burst). 
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6. Results 

The percentage of dust obscured afterglows for each model is computed summing all the 
absorption intervening in the four regions. Critical limiting magnitudes, beyond which 
the transient is no more visible, are taken from columns (5) and (6) of Tab. |21 while 
extinction curves used for standard Galactic dust and larger grain dust are respectively 
the dotted and the solid ones represented in Fig.^ For SMC extinction curve we use the 



one reported by Weingartner & Draine (2001) 



If in 5R2 there are only grains larger than 2^m (both in case of standard and graphite 
only grains) we use dashed curve of Fig. ^ but multiplied respectively by a factor Fi ~ 

a^^-^da/ /p^gpg a~'^'^da and F2 — 0.5 * Fi taking into account that total dust mass is 
decreased. The 0.5 factor takes into account that half of the dust population (silicates) 
has completely disappeared. 

In Tab. [SI results for different models are reported. It is evident that content in 3fi2 is 
almost not significant, whereas most of the extinction take place in K3, so the content of 
this region has a determinant role. In Fig. we plot percentages for relevant cases. 



12 Vergani, Molinari, Zerbi, Chincarini: Dust and dark GRBs 

Table 5. Optical (R) and infrared (K) transient loss (%) in the GRB-GMC association 
scenario. The models are described in text and Tab. 21 
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In the lower panels of Fig. [7| we show magnitudes distributions (R and K) of the 
simulated afterglows after 100s for GRBs at z—1. Together with the complete population, 
the distributions of observed transient for the different models are plotted. 

z=l 

80 1 \ \ \ \ 




dust model 

Fig. 6. Trends for the most relevant models of the percentage for prompt (RP, KP) and 
late (RT, KT) observations in R and K bands of dark bursts caused by dust extinction. 
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Fig. 7. Distribution of the 6 jet opening angles of our 5000 GRBs simulated sample 
(upper panel). The parent distribution follows a 0-^-^^ law after a fit to observed jet 
angles. The middle panel shows the translation of jet angles in R magnitudes (heavy line). 
The transformation lays on the zero point of a GRB990123-like event at z=l {9 = 0.086, 
R=10, K=7.5) and the assumption of the conservation of the total flux on the emitting 
solid angle (oc 6"^). Thinner lines are the luminosity function of the promptly observed 
(non absorbed in our simulations) afterglows in the R band. The four histograms refer 
to different models: 1, 3, Ismc and 3smc. In the lower panel the same set of distributions 
are plotted for the K filter (note that model 3smc does not differ from the original 
distribution, plotted as a heavy line). 

A comparison of our simulation thought in view of Swift and REM with observa- 
tional data available now would be misleading. From our point of view it is practically 
impossible to extract a statistically significant number (e.g. % of real dark bursts) from a 
set of observations performed with different telescopes (i.e. limiting magnitudes) and at 
different time from the burst. To stress this point we compute at different times the per- 
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centage of events for which the OT would not be observed even without dust extinction, 
considering the hght curves of our simulated GRBs and fixing R limiting magnitude of 
Riim = 20.5. Solid bold line in Fig. |S1 shows that for the majority of data observational 
times the percentage of lost events is significant. Dashed bold line of Fig. |H1 reports the 
computation made for the K infrared transient with A';,j^ = 19. 



Recently |Klose et al. (2003)| have dealt with a GRB (GRB020819) whose afterglow 
has been searched without success in K', 0.37 days after the burst down to K'=19 and in 
R band 0.13 days after the burst down to R=20.5. We report these observational times 
and magnitudes in Fig. |H1 we notice that we can not exclude that GRB020819 could not 
be a real obscured or high redshift burst. It is also worth noting that our computations 
remark that if the observations reach Rum — 24, the value at 0.13 days of the R band 
curve not affected by dust is null and that late K observations (more than 6 hours after 
the burst) are always sensibly biased. 

We add to Fig.|Hlthe R band results of our models (i.e. considering extinction) at the 
same magnitude limit. 



7. Conclusions 

We have considered different kinds of dust and their role in obscuring optical and NIR 
afterglows of GRBs simulated with different distributions inside a galaxy at z=l. In 
the case of GRBs occurring inside giant molecular clouds, our simulations show that 
dust in GRBs host galaxies has not the same properties of galactic dust, otherwise dark 
GRBs would be more than observational data say. Furthermore our results agree with 
the hypothesis of dust extinction curve of GRBs host galaxies similar to that of SMC. 
Moreover, if the host molecular cloud is characterized by large dust grains, high redshift 
plays a minimal role in the causes of dark GRBs. The opposite situation takes place if 
also inside the host molecular cloud dust has SMC properties. 

Once we will gather enough prompt K observations we will be able to estimate the 
nature of dust present in host molecular clouds from the fraction of lost K afterglows. 

We underline the fact that with present observational data it is impossible to produce 
a real statistic for dark bursts due to the late and varying time of observations and to 
the different magnitudes reached. 

In the future, from the afterglow results given by Swift and REM, we will be able to 
determine the nature of dark GRBs. Indeed with a further refinement of the model and 
a good statistics of bursts observed in various colors we will certainly be able to know 
which kind of dust is present in the environment of the burst. 
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t [sec] 

Fig. 8. Percentage of lost afterglows in function of observing time for GRB at z=l with 
our simulated light curves. Solid bold line represents the R band curve without consid- 
ering dust extinction, Rum, = 21; dashed bold line represents the K band curve without 
considering dust extinction, Kum = 19, whereas the other curves report R band results 
from our models calculated for Rum — 21. Two line intersections show the percentages 
in R and K relatives to the data of GRB020819. 
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